Spectral Density Scaling of Fluctuating Interfaces 
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Covariance matrices of heights measured relative to the average height of growing self-affine sur- 
faces in the steady state are investigated in the framework of random matrix theory. We show that 
the spectral density of the covariance matrix scales as p(A) ~ \~ v deviating from the prediction of 
random matrix theory and has a scaling form, p(X,L) = X~" f(X/L^) for the lateral system size L, 
where the scaling function f(x) approaches a constant for A -C and zero fori/ < A < X max . The 
values of exponents obtained by numerical simulations are v w 1.67 and cf> « 1.53 for the Edward- 
Wilkinson class and v w 1.59 and <f> m 1.75 for the Kardar-Parisi-Zhang class, respectively. The dis- 
tribution of the largest eigenvalues follows a scaling form as p(X ma x, L) = 1 / ' L h / ma i((A mal - L a ) / 'L b ) , 
which is different from the Tracy- Widom distribution of random matrix theory while the exponents 
a and b are given by the same values for the two different classes. 

PACS numbers: 05.40.-a, 02.50.-r, 68.35.Ct, 02.10.Yn 



Over recent decades, growth phenomena of fluctuating 
interfaces persists as a fascinating subject of statistical 
physics. Fluctuating interfaces are among the most well 
studied non-equilibrium systems due to their simplicity 
as well as ubiquity in nature and fundamental science 
[l| . Growth of interfaces governed by local rules typically 
lead to the formation of self-affine surfaces with universal 
scaling exponents for the surface width W(L, t) which is 
defined as the standard deviation of interface height over 
a system size L. The surface width characterizes the 
roughness of the interface and follows a scaling behavior 
W(L,t) = L a f(t/L z ), where the scaling function f(x) 
approaches to a constant for x 3> 1, and f(x) ~ x? for 
x <C 1 with the dynamic exponent z — a/ (3 2]. The 
exponents a, and /3 which are called the roughness and 
the growth exponent, respectively, which determine the 
universality classes of various fluctuating interfaces. The 
well-known universality class of the growing interfaces is 
the Kardar-Parisi-Zhang (KPZ) one which is predicted 
by the nonlinear Langevin equation [3j due to the slope 
dependent growth with the values of exponents a = 1/2 
and (3 — 1/3 for one dimension and was widely confirmed 
in numerical models [l[ ■ While lacking of nonlinearity in 
the growth process results in another universality class 
called the Edward- Wilkinson (EW) class where the val- 
ues of exponents are given by a = 1/2 and j3 = 1/4 for 
one dimension Q. 

The recent studies of fluctuating interfaces have dealt 
with other important characteristics beyond those for the 
scaling properties of the surface width. These include the 
distribution of the surface width and maximal and 
minimal height distributions [1,01, etc. Especially, the 
asymptotic distribution of KPZ height fluctuations for 



curved initial conditions has been computed exactly for 
solvable models [8| and found that it follows the Tracy- 
Widom distribution (TWD) [9] of a Gaussian unitary 
ensemble, which has been confirmed by a recent experi- 
ment on the electro-convection 1101 and the simulations 
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While for the flat initial condition, it is confirmed 
analytically [l2j and numerically [l3| that the KPZ height 
distribution exhibits the TWD of a Gaussian orthogonal 
ensemble. Although TWD has been first obtained in the 
statistics of the largest eigenvalues of random matrices 
belonging to the Gaussian ensembles @, it has been ap- 
plied to other areas [lij and its application to the growth 
phenomena has given the connection between the growth 
problem and the random matrix theory (RMT) [151 ] . 

RMT was initially proposed to explain the statistical 
properties of nuclear spectra [HI but the usefulness of it 
in understanding the statistical properties of a complex 
system makes RMT be applied in the various systems 
[lol \v?\ . The statistical properties of matrices with in- 
dependent random elements have been well described by 
the RMT and one can understand the statistical prop- 
erties of a system by comparing the spectral statistics 
of the system with the results of RMT. The empirical 
cross-correlation matrices appearing in the study of vari- 
ous complex systems such as the price fluctuations in the 
stock market [l8j], electro encephalogram(EEG) data of 
brain (l9| . variation of various atmospheric parameters 
[20l | . biophysical issues [2lJ and complex network [22j], 
have been analyzed in the framework of RMT. 

The random Wishart matrix that is one of the standard 
tools in the RMT is defined via the product W = j^XX^ 
of a M x N random matrix X having its elements drawn 
independently from a Gaussian distribution [23| . If X is 
a matrix whose elements represent some empirical data, 
then the Wishart matrix represents a empirical covari- 
ance matrix of the data, and the nondiagonal elements 
Wij of the covariance matrix have a direct interpretation 
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as cross-correlation coefficients between data AQ and Xj . 
Therefore if a certain complex system shows the spec- 
tral statistics same as those of the Wishart matrix, one 
may think there are not significant correlations, while if 
it shows the different properties from that, it could be 
regarded as there are some correlations. However, the 
random and the correlated properties are too entangled 
in a real complex systems to be simple elucidating the 
properties of correlations from the RMT analysis of cor- 
responding empirical covariance matrices. It informs just 
whether it is totally random or not. If the RMT is ap- 
plied to the problem of fluctuating interfaces in which the 
strong correlation of the variables of heights has been 
well understood would provide important insights into 
the universality of RMT for such a correlated system as 
well as novel criteria for statistical properties of a fluctu- 
ating interface. 

In this perspective, we would investigate further statis- 
tical properties of fluctuating interfaces under the frame- 
work of the RMT in this study. We construct the Wishart 
matrices by the product of the matrices having their ele- 
ments of relative heights obtained from two models which 
belong to the KPZ and the EW universality classes, 
respectively. The distributions of eigenvalues and the 
largest eigenvalues for each Wishart matrix are measured 
and compared with the results of the RMT. The ob- 
tained results have shown the totally different properties 
of those of RMT and we have found new scaling features 
of the distributions. 

Once the actual height hi(t) at the site i and the time t 
is generated, we define the relative height, Hu = hi(t) — 
(hi(t)), where the spatially averaged height (hi(t)) keeps 
on growing with time and by subtracting it, the relative 
height has the zero mean and the distribution of relative 
height reaches a stationary state in the late-time regime 
in a finite system. We consider the LxT height matrix H 
with elements Ha , where L is the lateral system size and 
T is the time interval we considered. We then compute 
the product symmetric matrix C — -^HH^ with elements 
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which represents the covariance matrix of heights and 
contains informations about height correlations between 
two sites. 

In the random matrix studies of eigenvalue spectra, 
the most popular property is the spectral density or the 
distribution of eigenvalues p(A). For the Wishart random 
matrix [23| . it was shown analytically that the spectral 
density p(A) is given by the Marcenko-Pastur(MP) law 
in the limit for M,N -> oo H: 
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FIG. 1: (a) The inset shows the spectral density of eigen- 
values p(A) of the covariance matrix of heights constructed 
by the EW model for L = 1024. The plot is the spectral 
density p(X, L) except for the largest eigenvalues for various 
L. (b) The distributions of the elements of the height ma- 
trices P(Hit,L) for the various L. The solid line represents 
the Gaussian curve, (c) The height correlations C(r,L) for 
various L. The solid line represents the curve of Eq. (d) 
The plot shows the data collapse of p(A, L) and the plot of 
(A 2 ) / (A) versus L is shown in the inset. 



where a is a standard deviation for elements of a M X N 
random matrix and m = M/N . The spectral density 
Pmp(A) shows the important features vanishing at both 
edges of the MP sea and exhibiting a sharp maximum 
near minimum edge. For finite M and N , the abrupt 
cut-off of p(A) is replaced by a rapidly-decaying edge |25j . 

To compare the spectral properties of the covariance 
matrix of heights with those of a random Wishart matrix, 
we have begun with the heights profiles generated by the 
EW model [4J which belongs to the EW universality class. 
We have obtained the numerical data for the one dimen- 
sional substrate with the system size L = 128, 256, 512, 
and 1024. Here we have only focused on the properties 
in the late-time regime in which the spacial correlation 
of heights is dominant and the correlation is independent 
on the time. Thus we counted the time interval T from 
the saturation time to the final time. 

The inset of Fig. [T] (a) shows the spectral density 
p(A, L) of the covariance matrix of heights constructed 
by the EW model for L = 1024. It shows the totally 
different feature from MP law. The spectral density p(A) 
follows a power-law behavior except for large eigenval- 
ues most which are due to the largest eigenvalues. So 
we considered the largest eigenvalue and the rest sep- 
arately. Figure [1] (a) shows the spectral density ex- 
cluded the largest eigenvalue for various system sizes. 
The straight guide line represents that the spectral den- 
sity obeys p(A) ~ \~ v with v = 1.67 ± 0.02. It is rea- 
sonable that the null hypothesis of no true correlation is 
rejected for the fluctuating interface which is a strongly 
correlated system. In addition to, it shows the difference 
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from the spectral density of the empirical covariance ma- 
trices most which have just appeared the distortion of 
the shapes of the bulk of spectral density [HI, [26| or ab- 
normal largest eigenvalues [l8j . The power-law behavior 
of the spectral density has been observed empirically in 
some systems such as the EEG data without the value of 
the exponent exactly measured (l9j. Analytical arg ument 
for it was proposed with the Levi matrices [28l [29| where 
matrix elements are distributed according to P(Xij) with 
P(Xij) ~ \Xij\-^ +l ^. For p > 2 the distribution has fi- 
nite variance while the variance diverges for < p < 2. 
In the case of the Wishart Levi matrices, the distribu- 
tion of eigenvalues has fat tails unlike the prediction of 
the MP law [13]. In Ref. [3l|, the multivariate Student 
distribution has been used as the power-like distribution 
and it has been obtained that the spectral density p(X) 
of the Wishart Levi matrices decays like \~^ l / 2+1 \ 

To check whether our spectral density follows the 
scheme of Levi matrix, we have measured the distri- 
butions of the elements of the height matrix P(Hn,L) 
for the various L. As shown in Fig. [1] (b) the dis- 
tributions P(Hu,L) has a scaling form, P(Hu,L) = 
1/L a g (Hit/L a ) and fall on a single curve, where g(x) 
is the Gaussian curve and a — 0.5. It results from that 
the only relevant scale is the roughness W being consis- 
tent with the scaling description in the late-time regime 
@,[32j]. Thus P(Hn) has no power-law tails, which means 
that the origin of the power-law behavior of p(X) is not 
in the fat tail of the distribution of matrix elements. And 
it shows that the MP law might be not valid any more 
even when the distribution of matrix obeys the Gaussian 
distribution. 

The elements of a Levi matrix are uncorrelated random 
variables, while in the late-time regime, the elements of 
the height matrix, i.e., the relative heights at different 
sites are strongly correlated as a following equation [32j 
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Figure Q] (c) shows the height correlations C(r, L) mea- 
sured by averaging the elements of the covariance matrix 
dj having r = \i — j\ which are excellent agreement with 
Eq. (HJ) for various system size L. It indicates that the co- 
variance matrix of heights has elements decreasing away 
from the main diagonal like Eq. (Q}. The long-ranged 
correlation of the elements might give rise to the power- 
law behavior of eigenvalue density p(X). It is comparable 
to the power-law random banded matrix (PRBM) model 
;33J which is defined as the ensemble of matrices with 
elements 
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where the matrix G runs over the GOE and a{r) ~ r~ a 
for large r. This exhibits an Anderson localization tran- 
sitions at a = 1 and allowed a detailed study of the 
wave function and energy-level statistics at criticality. 
Although the PRBM model represents ensembles of ma- 
trices with long-ranged off-diagonal disorder, its spectral 



properties are different from those of the covariance ma- 
trix of heights and thus it may serve as a another new 
critical random-matrix ensemble with long-ranged off- 
diagonal random hopping like Eq. 




-(b) „*» 



L=128 
L=256 
L=512 
L=1024 



10 o 

10 10 10 io 4 10 5 10' 



w 4 

" A 


o L=128 
□ L=256 
« L=512 
* L=1024 


; ;\ 

1 - 1 ^ta^r 





■0.05 0.05 0.1 
X 



FIG. 2: (a) The average of the largest eigenvalues {X m ax) 
as a function of the system size L for the EW model, (b) 
The distribution of the largest eigenvalues p(Xmax). (c) The 
standard deviation a(L) of \ ma x- (d) The data collapse of 
p(Xmax, L) for the various L. 

On the other hand, we found that the finite-size distri- 
bution of eigenvalues p(A, L) for the various system size 
L obeys a scaling form of the type, 



p(\,L) = \- v f 
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where A C (L) is a characteristic eigenvalue which scales 
as X C (L) ~ LP . f(x) is a scaling function satisfying the 
following properties: 



f(x) = const, for A < 

f(x) -> 0, for « A < A, ; 



(7) 



Figure [T] (d) shows the data collapse of the spectral den- 
sity p(X, L) with v = 1.67 and <fr = 1.53 and they fall on 
a single curve very well except for small eigenvalues. 

The exponent <f> can be measured by an alternative 
way. The n-th moment of the eigenvalue is obtained by 
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where A2 is the second largest eigenvalues and xi = 
X%jL^ . The integral has a finite value and we obtain 



(A") ~ ^{n+i-v) ^ 
Thus we can provide that 
(A 2 ) 
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We measured <fi = 1.53 ± 0.03 from it as shown in the 
inset of Fig. [T] (d) which is good agreement with the 
value used in the collapse of p(X,L). 

Next, we compared the properties of the largest eigen- 
values of the covariance matrix of heights with the results 
of the RMT. The MP law tells that the average of the 
largest eigenvalue (X max ) depends on the matrix size M 
like as (X max ) ~ M for large M and the typical fluctu- 
ations of X max are known to be described by the TWD 
9, 34]. We measured the average of largest eigenvalue 
^max(L) as a function of L and found that it scales as 
{^max) ~ L a with a « 2.09 (Fig. [2] (a)), which is also 
different from the result of the RMT. Also the distribu- 
tions of the largest eigenvalue p(\ mo .x) deviates from the 
TW curve as shown in the Fig. [5] (b) . Hence it does not 
follow the prediction of the RMT. 

The fluctuations of X m ax from its the average value 
come to be larger as L increases. We measured the stan- 
dard deviation u(L) of X max and found that cr(L) scales 
as cr(L) ~ L b with b ss 2.04. We rescaled the variable 
Xmax 'as x — (X max — L a ) / L b and obtained a scaling form 
of p(X max ,L) as follows, 
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FIG. 4: (a) The average of the largest eigenvalues {Xmax) 
as a function of the system size L for the RSOS model, (b) 
The standard deviation a(L) of Xmax- (c) The distribution 
of the largest eigenvalues p(X m ax)- (d) The data collapse of 
p(X m ax,L) for various L. 
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Figure [2] (d) shows that the distributions of largest eigen- 
values p{Xmax, L) fall on a single curve with a = 2.09 and 
b = 2.09 for the different system size L. Thus, the analy- 
sis of the fluctuating interfaces belonging to the EW uni- 
versality class under the framework of the RMT shows 
the results different from the conventional behaviors of 
the RMT and the new spectral scaling properties for the 
EW class. 
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FIG. 3: (a) The spectral density p(X, L) of the covariance 
matrix of heights constructed by the RSOS model except for 
the largest eigenvalue. The inset shows p(A) including the 
largest eigenvalue for L = 1024. (b) The data collapse of 
p{X,L) for various L. The inset shows the plot of (A 2 ) / (A) 
versus L. 



By applying the RMT to the analysis of the another 
universality class of the growth phenomena, the KPZ 
class, we would like to compare the spectral proper- 
ties between two universality classes. We constructed 
the covariance matrix of the heights generated by the 
restricted-solid-on-solid (RSOS) model j35[ which be- 
longs to the KPZ universality class. The spectral den- 
sity of this covariance matrix p(X) also follows a power- 
law behavior like the case of the EW model (Fig. [3] 



(a) ). However the value of exponent v was obtained as 
v = 1.59 ± 0.01, which is different from that of the EW 
model. The height correlation of the RSOS model also 
obeys the Eq. (j4j which indicates that the value of v 
depends on the details besides the height correlations. 
The spectral density p(X, L) for various L also follows 
the scaling form of Eq. (|6]). The data collapse of the 
spectral density p(X, L) for the RSOS model is shown in 
the Fig. [3] (b). The values of the exponent v and <f> were 
taken by 1.59 and 1.75, respectively. The inset of Fig. [3] 

(b) shows the plot of (A 2 ) / (A) versus L and we obtained 
(f> = 1.75 ± 0.03 of which value is also different from that 
of the EW model. 

The properties of the largest eigenvalues for the RSOS 
model are shown to be similar to those for the EW 
model. The average of the largest eigenvalues scales as 
(Xmax) ~ L a with a s» 2.08 (Fig. H](a)) and the standard 
deviation <j(L) of X max scales as a(L) - L b with b « 2.08 
(Fig. H] (b)). Figure |4] (c) shows the distribution of the 
largest eigenvalues which does not follow the TW curve 
like the EW model for the various system size L. The 
data collapse of p(X max , L) by the scaling form of Eq. 
(fTTjl falls on a single curve with a = 2.08 and b — 2.08 
(Fig. 11(d)). 

In summary, we investigated the spectral properties 
of covariance matrices of relative heights of fluctuating 
interfaces in the late-time regime. The spectral density 
of the covariance matrices follows the power-law behav- 
ior except for the largest eigenvalue, which is different 
from the MP law of RMT. It indicates that the random 
variables correlated like the relative height of a fluctuat- 
ing interface would give rise to the power-law behavior 
of the spectral density of the corresponding covariance 
matrix. As a finite-size effect the spectral density falls 
zero beyond the characteristic eigenvalue depending on 
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the lateral system size and has the scaling form. The 
values of exponents v and <j> related to the scaling form 
were given by the different values for the EW model and 
the RSOS model. It indicates that the statistical char- 
acteristics of the different classes of two models in the 
late-time regime are reflected to the spccral properties 
of them. The distribution of the largest eigenvalues also 
showed the different features from the TWD. The aver- 
age and the standard deviation of the largest eigenvalues 
follow the power-law behaviors with the lateral system 



size. The values of the scaling exponents do not give the 
difference between the EW model and the RSOS model. 
It would be desirable if it is applied to the behavior in the 
short-time regime and various interface models belonging 
to the different universality classes and the local spectral 
properties of the critical random matrix ensembles using 
the correlated heights are further studied. 

This work was supported by the 2012 research grant 
from Korea National University of Education. 
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